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Abstract 

Hartree-Fock calculations with the SKX Skyrme interaction are carried out to obtain alpha-decay 
Q values for deformed nuclei above 208 Pb assuming axial symmetry. The results for even-even nuclei 
are compared with experiment and with previous calculations. Predictions are made for alpha- 
decay Q values and half-lives of even-even super-heavy nuclei. The results are also compared for 
the recently discovered odd-even chain starting at Z = 112 and N = 165. 
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I. INTRODUCTION 


The existence and decay properties of super-heavy nuclei are one of the most fundamental 
problems in nuclear physics IS- There are now new data which confirm the existence of 
Z = 111 and 112 and their connection to lighter decay chains 0. The first data for Z = 114 
and Z = 116 also exist [[|], with suggested A = 288 and A = 292, respectively, but the 
A values are not certain since the connection to lighter nuclei is not known. Theoretical 
models for super-heavy nuclei have evolved from the macroscopic-microscopic models such 
as the finite-range droplet model with shell corrections ||, to fully microscopic deformed 
Hartree-Fock (HF) models such as those presented in [|6| and [7]. In addition to their intrinsic 
many-body nuclear structure importance, theoretical models for the prediction of the decay 
properties of the super-heavy nuclei are important when designing experiments since the 
techniques used will depend on the half-life and decay mode. 

In this paper we present a new set of Hartree-Fock results for alpha-decay Q-values of 
super-heavy nuclei. A global formula is used to calculate the half-lives. Our calculations 
are based upon a new computational program for solving the axial-symmetric HF equations 
and the new Skyrme interaction SKX ||. The reason for exploring results with another 
Skryme interaction beyond those used previously in 0 and 0 is that the alpha decay 
Q-value systematics are sensitive to the spherical and deformed shell-effects which depend 
upon the underlying parameters of the hamiltonians. There are several modern Skyrme 
parameter sets available, each of them determined with a different weighting and emphasis 
on the existing nuclear structure properties. 

The SLy4 parameters |J used in |J and the MSk7 0 parameters take into account overall 
spacing of the single-particle states such as those in 208 Pb. (In particular, it is common to 
constrain the Skyrme parameters to give an effective mass of unity which is required by the 
observed level spacing.) However, the SKX interaction || explicitly incorporates most of 
the observed single-particle levels in 208 Pb into the data set which was used to determine 
the parameter values. The single-particle energies for some proton and neutron particle 
states above the fermi surface in 208 Pb are compared in Table I. All of the calculations have 
some disagreement with experiment, however, SKX has the best overall agreement with 
experiment. It is important to explore the model-dependences of the HF results for the 
binding-energies of nuclei above 208 Pb. 
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In Sec. 0 we will briefly discuss the computation method for the binding energies and 
single-particle spectra for axial symmetric nuclei in the HF approximation. In Sec. H the 
results for alpha decay Q-values and lifetimes will be presented and discussed. The region 
of the chart of nuclides covered by our calculations is shown in Fig. [I]. This includes the 
region just above 208 Pb where the alpha-decay Q-values are measured, and extends up to the 
assumed spherical magic numbers of Z = 126 and N = 184 for the super-heavy predictions. 


II. METHOD OF CALCULATION 


There are several methods available for solving the Schrodinger equation for the single 
particle wave functions 4>i in a Skyrme-type potential with axial symmetry. Often a cylindri¬ 
cal coordinate system is chosen and the wave functions are expanded in a deformed harmonic 
oscillator basis with carefully adjusted oscillator strengths. Then the Hamiltonian is diago¬ 
nalized in an appropriately truncated space of basis states. Here, we solve the Schrodinger 
equation for fa in coordinate space with a spherical basis for the angular part of the wave 
function. In this approach there is a smooth transition to the case of spherical nuclei where 
the wave function simplifies considerably. 

The single-particle wave functions of protons (q = +1) and neutrons (q = —1) in a de¬ 
formed nucleus with axial symmetry are specified by three quantum numbers: the parity 
7r = ±1, the principal quantum number n = 1 , 2 ,..., and the projection of the total angu¬ 
lar momentum on the symmetry-axis = ±|,±|,.... We expand the wave functions in 
coordinate space 


(pq-irnflO ) _ ^ ^ ./'tj7rrmO(d')T 

K, 

with the radial wave functions f qn nKnfa) and vector-spherical harmonics 


( 1 ) 


y K n(r) = ^(/ m s v \ j ti)Y lm {f)x S v ( 2 ) 

mis 

which are obtained by coupling the orbital angular momentum l of the spherical harmonics 
Y] m and the spin s — \ of the spinors \sv to the total angular momentum j. The index 
n of the vector-spherical harmonic of Eq. (^|) specifies j = |/c| — ^ with l = hi — 1 for 
k > 0 and l = —k for n < 0. The sum in Eq. ([[]) runs over all hi with | k, |>| fl | +| where 
n G {1, —2, 3, —4,...} for positive parity states and hi e { — 1, 2, —3,4,...} for negative parity 
states. The Schrodinger equation for (p q7Tn n leads to a set of coupled differential equation 


3 


for the radial wave functions / g7rnre n(^) for all allowed k. In the actual calculation only 
contributions with |D| + | < \k\ < |Q| + ^ are considered. The radial wave functions are 
discretized on a grid with a step size of h = 0.2 fm inside an interval [0, R] with maximum 
radius R = (1.25 A 1 ^ + 12) fm for a nucleus with A nucleons. Derivatives are represented 
by five-point formulas. 

Particle densities g q , kinetic densities r q and spin-current densities J q appearing in the 
Skyrme-Hartree-Fock potentials and the energy density are easily calculated from the single 
particle wave functions 4> qnn o,. E.g., proton and neutron single particle densities are given 
by the multipole expansion 

Qq(r) = ^2QgL(r)Y L0 (f) (3) 

L 


where 


6qL{r) = 

7T nfi 


'MqTrntl 


r 


2 


^ kk' JqnnKflJqnriK'Q 

KK 1 


(4) 


with coefficients 


s~iLQ. 


dn yl n Y LO y K , n 


(5) 


The occupation probabilities w q7Tn a. in each state are determined by the BCS calculation. 
Only even values of L appear in the sum of Eq. (3) and contributions 0 < L < 10 are con¬ 
sidered in the calculation. Non-radial contributions of the spin-current density are neglected 
in the calculation. 

The binding energy of nucleus with A nucleons and Z protons in its ground state is 
calculated from 


BE (A, Z ) — — ( E m f + E pair — E cm — E rot ) 


( 6 ) 


with the mean-held contribution 


B-inf 



(7) 


which is obtained by integrating the Skyrme-Hartree-Fock energy density H{f) over the 
spatial coordinates. The pairing energy in the BCS approach is given by 

2 

( 8 ) 


Epair ^ ^ j | ^ ^ \J^■qTvni 1 (1 

q \nnfl , 


The pairing strength is G + 1 = 1.9 /\f~A MeV for protons and G-\ = 1.2/a J~A MeV for neu¬ 
trons, respectively. These values were obtained from a fit to experimental pairing gaps of 
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N = 146 isotones and Z — 92 isotopes. Only bound states are considered in the determi¬ 
nation of the occupation probabilities w q7rn n in the BCS calculation. The correction for the 
center-of-mass motion 

E cm = ^ (45 A~* - 25 A~^ MeV (9) 

is the same harmonic oscillator approximation as for spherical nuclei in the SKX parametriza- 
tion. The rotational correction is approximated by 

^rot. r,-j- 

Z,-L x 

where ^ is the many-body wave function of the nucleus in the BCS ground state. The 
moment of inertia I x for the rotation around the cc-axis is calculated in the cranking model. 



III. RESULTS 

The deformed HF calculations give the binding energies for the nuclei above 208 Pb. From 
these we calculate the alpha-decay Q-value, Q a = BE(A — 4, Z — 2) + BE( 4, 2) — BE(A, Z ), 
where the experimental value of BE( 4, 2) = 28.30 MeV is used for the alpha particle. The 
results for alpha-decay Q-values for even-even nuclei are shown in Fig. [2j. The purpose is to 
compare with the measured Q-values as well as to compare with the results based on the 
SLy4 interaction shown in Fig. 1 of Ref. ||. Comparison of the two figures shows that the 
results for SLy4 and SKX are remarkably similar even though they are based upon Skyrrne 
parameter sets which are determined completely independently. Both show good overall 
agreement with experimental Q-values [|1UJ to within a rms deviation of a few hundred keV, 
with the exception of a dip in the theoretical Q-values around N = 152 which is not present 
in experiment. The largest deviation for SKX is for the Q-value for 256 No at N = 154. 
We also show in Fig. |2| the comparison with experimental Q a values from the suggested 
placement of the Z = 116 — 114 — 112 decay chain Q. These also agree well with theory. 

Further results are shown in Fig. |3| as a function of neutron number and in Fig. | as a 
function of proton number. The points are connected in these figures for a given N — Z 
value in order to emphasize how the Q-values changes in a given decay chain. Comparisons 
are made to the finite-range droplet model (FRDM) || and to the deformed HF-BCS calcu¬ 
lations based on the MSk7 Skyrrne interaction 0. The results for SKX and MSk7 are very 
similar even for the extrapolation to large N and Z. The FRDM results are similar to the 
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HF in the region where data are available but become more different for the extrapolation 
to heavier nuclei. 

Much of the data for the super-heavy nuclei are for odd-even decay chains. These are 
more difficult to calculate and compare with experiment since the deformed level density is 
high and the observed nuclei may be in isomeric states. These must be considered carefully. 
For this paper we compare with the Q-values observed for the recently confirmed decay 
chain for N — Z = 53 starting at 2 ' 7 112 in Fig. || In the calculation we assume that the 
nucleus is in its lowest energy deformed single-particle state. The results are also compared 
to the FRDM and MSk7 models. As in Figs. [| and [|. the SKX and MSk7 results are close 
to each other and both are close to experiment, with perhaps SKX being in best agreement 
with experiment. The FRDM results do not agree as well in detail with experiment. In the 
deformed HF the jump in Q-value between N — 161 and N = 163 observed in Fig. |5] comes 
from a deformed shell gap at N = 162 and Z = 108. These deformed gaps are also found 
with the SLy4 interaction ||. We note the semi-empirical shell-model mass approach [11 


cannot and does not take into account these deformed shell gaps and cannot reproduce any 
of the fine structure in the Q-value systematics. 

The distribution of the single-particle energies for protons and neutrons is shown in Fig. 
as a function of the neutron number N in even-even nuclei for the N — Z = 60 decay chain. 
Nuclei with small N are well-deformed and become more and more spherical with increasing 
N. Proton and neutron shell gaps are readily seen in both spherical and deformed nuclei. 
The proton Fermi energy increases smoothly with increasing N and becomes positive for the 
nucleus with N = 188 which is predicted to be proton-unbound in the SKX parametrization. 
The neutron Fermi energy decreases only slightly with increasing N. It crosses the well 
defined shell gap at N = 162 which was mentioned above. 

The alpha-decay half-life is important for determining how the alpha decay of super-heavy 
nuclei competes with fission. The extrapolated half-lives are also important for choosing 
the type of experimental techniques used for their identification. The main theoretical 
uncertainty for the calculation of the assumed L — 0 decays of even-even nuclei is in the 


alpha-decay Q-value. To calculate the half-lives we use the empirical result obtained in [12 


log 10 [T 1/2 /s] = 9.54 (Z - 2) 0 ' 6 /v 7 Qa/MeV - 51.37 


( 11 ) 


In Fig. we show the half-life calculated from Eq. ([IT]) and from the experimental Q-values 
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T0|1 • The results are compared to the experimental half-lives. The excellent agreement 


between experiment and theory shows that preformation and decay systematics implied by 
Eq. are adequate for a determination of the alpha-decay half-life to within about a 
factor of three. 

The predictions for the half-lives of heavier nuclei based upon the theoretical Q a values 
from our SKX calculations results are shown in Fig. ^|. The agreement with experiment is 
satisfactory except with near N = 152 where the kink in the experimental half-lives is not 
reproduced by the theory. The experimental half-lives for the suggested placement of the 
Z = 116 decay chain [Q are also in reasonable agreement with theory. One observes an 
island of relative stability starting at IV = 164 where the half-lives for Z & 53 remain at the 
msec level or longer until N & 174 where they start to become shorter. 


IV. SUMMARY 

We have presented a new calculation for the alpha decay Q values for super-heavy nu¬ 
clei based upon deformed Hartree-Fock calculations with the SKX Skyrme interaction. A 
new computational method is used to carry out axially-symmetric deformed calculations. 
Agreement with experimental data including the recently observed Z = 112, Z = 114 and 
Z = 116 decays is obtained at the rms level of a few hundred keV. Deformed shell gaps at 
N = 162 and Z = 108 lead to jumps in the Q values which are consistent with experiment. 
The Q values have been used to calculate alpha-decay half-lives which are in reasonable 
agreement with theory. Predictions for the Q values and half-lives up to the proton drip 
line at N = 184 and Z = 126 are made. 
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TABLE I: Single-particle energies for states in 208 Pb and the rms difference between experiment 
and theory. 


single-particle energy (MeV) 

orbit exp SKX SLy4 

MSk7 

Trlhg/2 —3.80 —4.26 —3.82 

-3.47 

7r2.7y2 —2.90 —2.83 —2.93 

-3.19 

7rl* 13/2 -2.29 -2.18 -1.46 

-2.62 

7r2/ 5 / 2 -0.98 -0.70 -0.37 

-0.93 

v2g 9/2 -3.94 -3.46 -3.14 

-4.05 

i/l*n /2 -3.16 -2.81 -1.57 

-2.02 

i/lji5/2 -2.16 -1.97 -0.64 

-2.53 

v2g 7/2 -1.45 -1.00 0.08 

-0.97 

rms 0.34 1.05 

0.51 
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FIG. 1: Chart of the nuclides for even-even nuclei with N > 138 and Z > 86. Solid squares denote 
nuclei with experimentally known masses. Nuclei discussed in this paper are located on the solid 
lines with constant N — Z (indicated by the numbers). Open circles indicate proton-unbound nuclei 
in the SKX parametrization. 
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FIG. 2: Q-value for a-decay as a function of the neutron number N connected by lines for the 
given Z values. The experimental values are shown by the solid circles. The results from the SKX 
deformed HF calculations are given by the open circles. 
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FIG. 3: Q-value for a-decay as a function of the neutron number N for even-even nuclei located on 
the solid lines of Fig. |]. Predictions from the finite range droplet model (open diamonds: FRDM), 
and two Skyrme Hartree-Fock parametrizations (open squares: MSk7 and open circles: SKX) are 
compared with the experimental data (solid circles). Different decay chains with the same value of 
N — Z (indicated by the numbers) are connected by solid lines and shifted vertically by the amount 
(in MeV) shown in parentheses. 
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FIG. 4: Same as Figure | as a function of Z. 
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FIG. 5: Same as Figure |3] but for the nuclei with N — Z = 53. 
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FIG. 6: Single particle energies of protons (a) and neutrons (b) above —20 MeV in even-even 
nuclei with N — Z = 60 in the Skyrrne Hartree-Fock calculation with the SKX parametrization. 
The Fermi energies of protons and neutrons are denoted by open circles. 
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FIG. 7: Half-lives of even-even nuclei as a function of the neutron number N. Open symbols 
indicate half-lives calculated with Eq. (4) with the experimentally measured Q a values or with Q a 
values from the Audi-Wapstra mass extrapolation. Solid symbols denote experimental half-lives. 
Decay chains with constant N — Z value are connected by solid lines. 




FIG. 8: Half-lives of even-even nuclei as a function of the neutron number N. Open (solid) symbols 
indicate half-lives calculated with Q a values from the SKX parametrization (experiment). Decay 
chains with constant N — Z value are connected by solid lines. The experimental half-lives from 
the suggested placement of the 292 116 decay chain are shown by the cross-filled boxes. 
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